Tunable Whispering Gallery Mode Resonators for Cavity Quantum Electrodynamics 
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We theoretically study the properties of highly prolate shaped dielectric microresonators. Such 
resonators sustain whispering gallery modes that exhibit two spatially well separated regions with 
enhanced field strength. The field per photon on the resonator surface is significantly higher than 
e.g. for equatorial whispering gallery modes in microsphere resonators with a comparable mode 
volume. At the same time, the frequency spacing of these modes is much more favorable, so that a 
tuning range of several free spectral ranges should be attainable. We discuss the possible application 
of such resonators for cavity quantum electrodynamics experiments with neutral atoms and reveal 
distinct advantages with respect to existing concepts. 

PACS numbers: 42.55.Sa, 03.65.Ge, 42.50.Pq 



Optical microresonators have attracted much interest 
in the field of linear and nonlinear optics as well as in 
cavity quantum electrodynamics, as summarized in the 
recent review by K. Vahala Q. One of the reasons for 
this considerable attention is the prime importance of a 
small mode volume V for such applications. So far, opti- 
cal microresonators based on Fabry-Perot cavities, pho- 
tonic crystals, and whispering gallery modes have been 
realized [l[. I n addition to V, certain applications, like 
e.g. cavity quantum electrodynamics Q, require also a 
minimization of the ratio V/Q, where the quality factor 
Q is determined by the energy storage time in units of 
the optical period. 

Very high Q-values have been realized for Fabry-Perot 
cavities with ultrahigh reflectance mirrors p|. However, 
the mirror components are involved and costly to manu- 
facture. Furthermore, due to the modular construction, 
such resonators require an elaborate active stabilization 
of their resonance frequency. Photonic crystal based res- 
onators have a much better passive stability and yield 
extremely small mode volumes. Even though the experi- 
mentally realized Q-values fall well below the theoretical 
limits, they therefore reach record values for V/Q fl. 
However, the problem of tuning the resonance frequency 
of such structures has not been solved to this date. 

The highest Q-values for optical resonators to date 
have been realized with whispering gallery modes 
(WGMs) in fused silica microspheres p, where con- 
tinuous total internal reflection confines the light to a 
thin equatorial ring near the surface of the sphere. In 
combination with their relatively small mode volume of 
V ~ 1000 A 3 , where A is the optical wavelength, these 
modes are therefore ideal candidates for cavity quantum 
electrodynamics experiments. 

Optical probing and output coupling requires phase 
matching of the WGMs to propagating light fields. This 
has been achieved with various methods, all of which use 
an auxiliary dielectric structures (prisms 6] , eroded or ta- 
pered optical fibers 0, 0| , or fiber /prism hybrid systems 
P) that are introduced into the external evanescent field 



of the WGM in the vicinity (< A) of the sphere surface. 
Using such a set-up, strong coupling between atoms in a 
dilute vapor and a microsphere WGM has been realized 
jldj l . However, in spite of considerable research activities, 
a controlled strong coupling of a suitable dipole emitter 
to a microsphere WGM has not yet been accomplished. 
The main difficulty to be overcome consists in control- 
lably placing the emitter as close as possible to the di- 
electric surface. In particular, due to the two-dimensional 
character of equatorial WGMs, the prism or fiber coupler 
severely limits the mechanical and optical access. 

Another practical difficulty stems from the limited tun- 
ability of microsphere WGMs. For a typical microsphere 
of « 50 /im, tuning over one free spectral range (FSR) 
requires a change of the optical path length of the WGM 
of about 4 x 10 -3 . On the other hand, the tempera- 
ture dependence of the refractive index of silica is only 
dn/dT « 1.3 x 10 -5 K _1 . This index change is the dom- 
inant temperature effect in silica. Therefore, a temper- 
ature variation can only be used for fine tuning in this 
case. The only practicable solution for tuning is thus to 
elastically deform the sphere through mechanical strain. 
For a ~ 80 /am microsphere, tuning over one half of 
an FSR has been achieved with this method, limited by 
the mechanical damage threshold of the resonator set-up 
[llj . In addition, for a w 50 jam microsphere with 
Q « 10 9 the resonance linewidth is only 10 7 of an FSR. 
Therefore, even if such a tuning device reached a full FSR 
tuning range, it would be a challenging task to guarantee 
the corresponding passive stability. 

Here, we propose the use of WGMs in highly pro- 
late dielectric resonators with a cylindrical symmetry. 
Of course, such a structure also sustains "equatorial" 
WGMs, i.e. with a corresponding ray path that is closed 
after one revolution around the resonator axis. In the 
following, however, we will consider modes for which the 
light spirals back and forth along the resonator axis be- 
tween two turning points, separated by 2z c , see fig. [I] 

From this very simple ray path picture it is already ap- 
parent that such modes exhibit two spatially separated 
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FIG. 1: Sketch of the resonator geometry. The ray path corre- 
sponding to the whispering gallery modes under consideration 
is schematically indicated. Note that the thickness variation 
of the resonator structure along the z-axis is amplified for 
illustration purposes. 



caustics located at ±z c with an enhanced field strength. 
Note that, in analogy to a charged particle in a magnetic 
bottle, the light is in fact confined along the z-axis be- 
tween the two turning points by an angular momentum 
barrier. For this reason, we have adopted the denomi- 
nation "bottle modes" from , where such a resonator 
shape was theoretically constructed from an equidistant 
spectrum of eigenmodes using WKB approximation. 

We will now present a full scale wave equation calcula- 
tion for such a resonator shape and theoretically explore 
important resonator properties including field distribu- 
tion, mode volume, and tunability. We consider a res- 
onator profile which is approximately parabolic along z: 



R(z) « Rq 1 
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(1) 



8 /im as the maximal radius and 



5 /im 2 for the curvature of the profile. We 



where we take Ro 
Ak 2 * 10 

use the method of adiabatic invariants well known 
from the description of the dynamics in magnetic bot- 
tles. Indeed, due to the small variation of the resonator 
radius (amplified in Fig. [I}, dR/dz <C 1, and since we are 
considering maximum angular momentum modes located 
close to the surface of the resonator, the radial compo- 
nent k p of the wave vector is negligible with respect to 
k z and k^. Therefore, the total wave number is 



27rn/A , 



(2) 



where n and A are the effective index of refraction 
and the wavelength in vacuum, respectively. Now, due 
to cylindrical symmetry, the projection of the angu- 
lar momentum onto the z-axis is a conserved quantity, 
d z k(p(z)R(z) = 0. Furthermore, the axial component of 
the wave vector vanishes in the caustics, k z (±z c ) = 0, so 
that k(p(^zz c ) = k. The azimuthal and axial components 
of the wave vector thus read: 



k<p(z) = kR c /R(z) and 



k z (z) = ±kJl-(R c /R(z)Y 



(3) 
(4) 



where — z c < z < z c and R c = R(z c ). 

In order to check that the paraxial approximation of a 
negligible radial wave vector component is well justified 
for the resonator geometry of eq. (QJ, we use the relation 
k p — (dR/dz)k z . For z c ~ 70 jam this yields \k p (z)\ < 
6xl()- 3 \k z (z)\ in the region — z c < z < z c . Inserting this 
value as a correction into eq. we estimate the error 
in eq. (@J) to be smaller than 2 x 10 -5 . 

Using the adiabatic approximation, we will now estab- 
lish the wave equation and determine the eigenfunctions 
of the resonator. Due to the cylindrical symmetry, the az- 
imuthal part of the wave equation can be separated with 
solutions proportional to exp(zra<^), where m is the az- 
imuthal quantum number. The solution can thus be writ- 
ten z) exp(irrnp). The radial quantum number p will 
be fixed to its minimum value p = 1, corresponding to 
modes located at the surface of the resonator. Now, using 
the adiabatic approximation along the z-axis, z) can 
be separated in a product of two functions <l>(p, z)Z(z), 
where & is the solution of a Bessel equation, 
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P dp 
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Here, k^ is given by eq. (0) and, using kR c — m, equals 
k(p(z) — m/R(z). Insertion into the above equation yields 



d 2 $ 



p dp 
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In order to determine the radial wave function from this 
equation, we first have to fix the boundary conditions at 
the resonator surface, which depend on the polarization 
of the mode. We consider a TM mode, i.e. linear polar- 
ization with the electric field parallel to the z-axis. In 
this case, the solution of equation © can be written 



A rn J rn (k i pp) 



(1) /fcypx 
V n / 



(p < R(z)) 
(p > R{z)) 



(7) 

(1 2) 

where J m and are the Bessel and Hankel functions, 

respectively. The coefficients A m and S m are determined 
such that they satisfy the matching conditions for <3> and 
its derivative, d^/dp. 

The differential equation for Z also depends on the res- 
onator profile R(z). In order to find an analytic solution, 
we choose the explicit form R(z) = Ro/ y/l + (Akz) 2 
which, given that (Akz) 2 < 0.05 in our case, realizes 
the parabolic profile of equation (QJ to a very good ap- 
proximation. Using this form, the problem reduces to a 
simple harmonic oscillator, 



d 2 Z 

dz 2 



i 2 m 2 Ak 2 



Rl 



Z = 0. 



(8) 



In analogy to the harmonic oscillator problem, we can 
therefore identify the total and the potential energy as 



3 




E = k 2 - m 2 /Rl and V(z) = (AE m z/2) 2 , respectively, 
where we have set AE m = 2mAk/Ro. Furthermore, the 
condition that Z be square integrable leads to a discrete 
set of energy levels E mq = (q + 1/2) AE m , where q is the 
axial quantum number, corresponding to the number of 
nodes of the wave function along z. This allows us to 
deduce the allowed eigenvalues 

k mq = [m 2 /Rl + (q+ 1/2) AE m ] 1/2 , (9) 
with the corresponding solution of equation |JBJl, 

Z mq {z) = C mq H q (^l^^e W [-^z 2 ) , (10) 

where H q is the Hermite polynomial of order q with the 

normalization constant C m q = [AE m / (7r2 2q + 1 (ql) 2 )] 1 ^ . 

From the wave function \I/ mg (p, z) = & m (p, z)Z mq (z) 
we can now calculate the intensity distribution 
Imq(p,z) oc \Sfr mq (p, z)\ 2 of the bottle modes. We con- 
sider the case where the wavelength X mq matches the 
D2 line of Cs at 852 nm. We choose the resonator ra- 
dius at the caustics to be R c = 7.8 //m, compatible 
with a quality factor of Q « 10 9 [l4|. Moreover, we set 
Ak = 0.0032 /im -1 , q = 80, and m = 76, i.e. the highest 
possible m value. Figure shows I mq (p, z) for the cor- 
responding bottle mode. The intensity in the caustics, 
located at z c = ±69 /im, is enhanced by a factor of four 
compared to the peak value around z = 0. 

The intensity distribution directly yields the mode vol- 
ume V of the bottle modes. It is obtained by normalizing 
Imq(p, z) to unity and by integrating it according to: 

Vmq = fff P0( 'e{ p , z )Il!i^l p dpd t pdz, (11) 

J J JO 1 mq 

where s(p, z) = n 2 inside the resonator and 1 outside, n 
being the index of refraction of the resonator material. 
The upper limit for the radial integral, po(z), is chosen 
such that it coincides with the first zero of the radial wave 
function outside the resonator, <£ m (po, z) = 0, in order to 
include the effects of the evanescent field. For the bottle 
mode of fig. we find V76,80 = 690 /im 3 . This value is 
only twice as large as the mode volume of a microsphere 
with « 50 pm. 



FIG. 3: Coupling strength of a Cs atom to the mode of the 
bottle resonator as a function of the distance from the res- 
onator surface (solid line). The atom is placed at the point 
of highest field intensity along the the resonator axis, z — z c . 
For comparison, the coupling to an equatorial WGM of a 
= 50 pm microsphere is also given (dashed line). Due to its 
smaller mode volume, the maximum coupling coefficient for 
the microsphere is a factor of 1.4 larger. However, the cou- 
pling at the resonator surface is 1.5 times larger for the bottle 
resonator, because its evanescent field is more pronounced. 



Now, the maximum coupling strength of a given dipole 
emitter to a single photon in the resonator mode is in- 
versely proportional to the square root of the mode vol- 
ume, g cx 1/y/V. However, for WGMs we can only access 
the field outside of the resonator. The actual parameter 
of interest is therefore the coupling strength at the sur- 
face of the resonator. We have calculated this value for 
our bottle resonator at the position of the caustic z c . For 
the D2 transition of a Cs atom we find g/2ir w 90 MHz. 
In spite of the larger mode volume, this exceeds the ex- 
pected coupling strength to the equatorial WGM at the 
surface of a ^50 pm microsphere by a factor 1.5. 
The physical reason for this result lies in the fact that, 
due to the smaller diameter of the bottle resonator, its 
evanescent field is more pronounced than for the larger 
microsphere (see fig. 0J). 

We now examine the mode spectrum of the bottle res- 
onator, given by eq. (0. For our geometry, i.e. RoAk <C 1 
and m ~ one finds that the frequency spacing between 
modes for which the azimuthal quantum number m dif- 
fers by one approximately equals Av m « c/2i\nR^ « 
4 THz, just as large as for a microsphere resonator of 
the same radius. However, the mode spacing concerning 
the axial quantum number q is much smaller and is given 
by Av q w cAk/27T w 150 GHz, one order of magnitude 
smaller than for a ^ 50 pm microsphere. Tuning over 
one free spectral range therefore only requires a change of 
the optical path length of 4 x 10 -4 . This can be realized 
by varying the resonator temperature by 30 K, a param- 
eter range that is easily accessible in an experiment. 

Mechanical strain could also be used to tune the bottle 
resonator. Concerning its mechanical properties, we ap- 
proximate the resonator as a cylinder of radius Rq. This 
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is well justified in our case and allows an easy estimation 
of the effect of the strain on the mode frequency v mq . 
Using the Poisson coefficient and the elasto-optic coeffi- 
cients of silica 15] and neglecting the modification of the 
curvature Ak along the resonator axis, we find 

Avrnq/vmq « -AR/R - An/n « -0.2 AL / L (12) 

for the TM mode (polarization parallel to the strain), 
where L is the length of the resonator. One FSR tuning in 
our case thus requires a length change of about 2 x 10 -3 . 
Using the value of 7.2 x 10 10 Pa for the elasticity module 
of silica, this implies a strain of about 0.15 GPa, more 
than one order of magnitude smaller than the typical 
damage threshold for silica fibers of 3 GPa Strain 
tuning over several FSR therefore seems possible. 

Apart from the advantageous properties concerning 
mode volume and tunability, the bottle resonator geom- 
etry and the resulting spatial structure of the modes also 
offer interesting advantages with respect to their han- 
dling and possible applications. Since bottle resonators 
exhibit two spatially well separated regions of enhanced 
field strength (see fig-EJ, they are true two-port devices, 
contrary to microspheres. One of these two caustics could 
therefore be used to couple light into and out of the res- 
onator by means of a thin unclad optical fiber |g, while 
a dipole emitter could be coupled to the other caustic. 

In practice, it has been demonstrated that a prolate 
resonator geometry with the dimensions considered here 
can be realized on unclad optical fibers usin g fib er pulling 
techniques and CO2 laser microstructuring 17]. In order 
to obtain the desired high Q values, these processes will 
have to be optimized to yield a surface quality that mini- 
mizes scattering losses. Furthermore, the bulk absorption 
of the fiber cladding, which would exclusively guide the 
light in this case, has to be sufficiently low. Using low- 
loss telecommunication fibers, we are however confident, 
that this condition can be fulfilled. 

Summarizing, we have studied the properties of whis- 
pering gallery modes in highly prolate shaped silica res- 
onators. We have shown that such "bottle resonators" 
sustain modes that exhibit two spatially well separated 
caustics with an enhanced field strength. For a distance 
of 140 /im between the caustics and a resonator diameter 
of 16 /im we find a mode volume of 690 /im, only twice 
as large as for an equatorial whispering gallery mode in 
a silica microsphere with a diameter of 50 /im. In spite 
of this larger mode volume, the coupling strength of an 
atom to the evanescent field of the mode outside of the 
bottle resonator is larger than for a microsphere: Due 
to the smaller radius of the bottle resonator, its evanes- 
cent field is stronger. For the D2 transition of Cs, we 
calculate a coupling strength of g/2ir ~ 90 MHz, much 
larger than the atomic line width and the expected cav- 
ity line width for an estimated quality factor of 10 8 -10 9 . 
The bottle resonator should therefore allow to enter the 
strong coupling regime in neutral atom cavity QED. 



At the same time, we have shown that the mode spac- 
ing for the bottle resonator is one order of magnitude 
smaller than for a microsphere of 50 /im diameter. By 
varying the temperature or applying mechanical strain 
to the resonator, tuning over more than one free spec- 
tral range should therefore be possible. Moreover, the 
fact that bottle resonators are true two-port devices is 
an important advantage when it comes to probing the 
light field in the resonator mode while simultaneously 
coupling a dipole emitter to it. Finally, it has recently 
been proposed to trap and guide atoms around thin un- 
clad optical fibers using a two-color evanescent light field 
[l^ . Our resonator design seems particularly suited to 
be combined with such a surface trap in order to couple 
trapped cold atoms to the resonator mode. 
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